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v^ Abstract. Using the Multipolar post-Minskowskian formal- 

\ I ism, we compute the linear momentum flux from black-hole 

^ binaries in circular orbits and having no spins. The total linear 

momentum flux contains various types of instantaneous (which 
are functions of the retarded time) and hereditary (which de- 
pends on the dynamics of the binary in the past) terms both 
of which are analytically computed. In addition to the inspiral 
contribution, we use a simple model of plunge to compute the 
kick or recoil accumulated during this phase. 



1. Introduction 

In classical radiation theory, any form of radiation that has a preferred 
direction (anisotropic emission) results in the recoil of the system in 
the opposite direction. Such a recoil has important consequences 
in astrophysics like the pulsar acceleration due to the radiation 
asymmetry [IT]. This effect can be understood in terms of multipole 
moments of the system, more specifically interference between different 

* Present address: USER Thiruvananthapuram, Computer Science Building, 
College of Engineering Campus, Thiruvananthapuram, 695016, India 



multipoles [17]. It was argued in Ref. [17], on general grounds, that 
the leading recoil effect in the case of electromagnetic radiation could 
come from the interference between electric dipole and electric quadrupole 
or magnetic dipole (See Ref. [12] for a more detailed discussion). A 
similar effect can happen in the case of gravitational radiation [17 . Since 
the leading emission in the case of Gravitational Waves (GWs) in GR 
is quadrupolar, the lowest order recoil effect in GR arises from the 
interference of mass quadrupole and mass octupole or current quadrupolej^ 
This GW recoil can be seen as a consequence of linear momentum flux 
emission by the radiating system. 

Let us now focus on a particular type of GW source known as coalescing 
compact binaries (CCBs). As the name indicates, these are systems 
consisting of neutron stars (NS) or black holes (BH) which go around each 
other in a bound orbit. The orbit keeps on shrinking due to gravitational 
radiation reaction until the two objects merge to form a compact object, 
most likely a BH. As is well known, the GW emission is not isotropic 
(due to the quadrupolar nature) and thus the emission is beamed in some 
direction. As the two bodies move in their respective orbits, this direction 
also keep changing. If the orbit is closed, then over an orbit, the binary's 
centre of mass would return to the same point where it started, making 
the net recoil zero. However, if the orbit is not closed, as is the case for 
CCBs, there can be an accumulation of recoil over these orbits until the 
two bodies coalesce into a single BH, eventually manifesting as a kick to 
this newly formed BH. This is the GW recoil in the context of CCBs. 

Gravitational recoil has important implications for astrophysics, 
especially models of black hole formation and growth (see e.g. p[F). If 
the recoil velocity of a merging compact binary system is greater than the 
escape velocity of the galaxy, it can result in the ejection of the newly 
formed black hole from the galaxy. Hence very accurate estimation of the 
recoil velocity of the merging compact binaries is important to understand 
the astrophysics of black hole formation and growth. 

Evolution of CCBs can be divided into three phases: adiabatic 
inspiral, nonlinear merger and perturbative ringdown. The inspiral 
part and ringdown phases can be accurately modelled analytically using 
post-Newtonian (PN) approximation to general relativity [3] and BH 
perturbation theory |19| , respectively. Due to recent successes in numerical 
relativity (see Refs. [THJ E] for reviews), the highly nonlinear merger 
phase can be modelled quite accurately by numerically solving Einstein's 
equations. 

In this work, we use the PN approximation to calculate the linear 

^ Both mass octupole and current quadrupole have the same parity. 



momentuni loss and resulting GW recoil from a black- hole binarjjj moving 
in circular orbits up to 2.5PN order extending the earlier works of Fittchett 
(Newtonian order) [TO], Wiseman (IPN order [5Tj) and Blanchet et al. 
(through 2PN [5]). A more detailed account of the contents of this paper 
can be found in Ref. |16) . 

This article is organized in the following way. In Sec. [2J we write 
down the multipolar expansion of the linear momentum flux and the final 
expression for the flux in terms of source multipole moments. Section [3] 
derives the 2.5PN expression for the recoil and Sec. |4] discusses the 
numerical estimates of recoil. Conclusions are given in Sec. [5] 

2. Multipole Expansion for Linear momentum flux 

The general formula for linear momentum flux in the far-zone of the source 
in terms of symmetric trace-free (STF) radiative multipole moments is 
given in ^2Qj and at relative 2.5PN order it takes the following form (see 
e.g. Eq. (4.20') of Ref. [20.) 
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In the equation above Ul and Vl denote the mass and current type 
radiative multipoles of the source (suffix L captures the multi-index 
structure of the U and V moments) and Uj^ and ¥£" denote the n"^ time 
derivatives of the moments. 0{l/c^) denotes the omission of terms higher 
than 3PN w.r.t the leading term. The moments which appear in the above 
formula are functions of retarded time T— — , where R denotes the distance 
of the source relative to the observer and T the time of observation, both 
in radiative coordinates. In the MPM formalism Ul and Vl are related 
to canonical moments Ml and Sl (Eqs. (5.4)-(5.8) of \A\) which in turn 
are related to source moments {II, Jl, Xl, Wl,Yl, Zl} (Eqs. (5.9)-(5.11) 
of [1]). Using these inputs, one can re-express the radiative multipole 

® Note that the expressions for the linear momentum flux (Eq. (|2|) and that 
for the recoil velocity (Eq. nj) are even applicable to compact binary systems 
involving NSs as components. However, the plunge computations presented here 
assume that both the components of the binary are BHs. 



moments in terms of source moments. Then one notices that the total 
hnear momentum flux consists of two parts: one type of terms are functions 
of the retarded time called instantaneous terms and the other which are 
sensitive to the dynamics of the source in its entire past called as hereditary 
terms. Explicit expressions for the two types of terms is given in Eqs. (2.3)- 
(2.5) of Ref. [16j . Further, the explicit expressions for various source 
multipole moments are given in Eqs. (3.1)-(3.4) of Ref. [TB]. The only 
other input we require is the 2.5PN accurate equations of motion which 
can be found in Ref. [HHJ]. Using these inputs, and working in harmonic 
coordinates, we obtain the total linear momentum flux in terms of gauge 
independent variable x ~ (muj)'^'^ as 
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In the equation above n and A are related to the phase angle tp as 



n = cos -0 Gx + sin -0 ey , 
A = —sin Gx + cos Gy , 



(3) 
(4) 



3. Computation of the RGCoil 

Having computed the linear momentum flux for compact binaries in 
circular orbits, one can obtain the recoil velocity by integrating the 
momentum balance equation 
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Performing the integration we find, 
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4. Numerical estimates of the recoil velocity including plunge 
contribution 

As is evident from Eq. ([7]) , the recoil velocity depends on the (symmetric) 
mass ratio v and does not depend on the total mass. This is consistent with 
our understanding that the origin of the recoil is in the mass asymmetry. 
We finally want to numerically estimate the recoil velocities as a function 
of V. Our calculation based on 2.5PN approximation yields a maximum 
kick velocity of ^ 4km/sec as opposed to ^ 22 km/sec of the 2PN model 
given in Ref. [5 . Thus 2.5PN estimates predict a smaller kick velocity 
than the 2PN model. This result is obtained by integrating the linear 
momentum flux till the Innermost Stable Circular Orbit (ISCO), up till 
which PN approximation is considered to be valid. 

Since the dominant contribution to the recoil comes towards the late 
stages of the binary evolution, we incorporate the contribution from the 
'plunge' phase of the evolution beyond the last stable orbit, following and 
extending a model proposed in Ref. [5] . The method may be considered to 
be less sophisticated version of the Effective One Body approach [6, 7 . In 
this model, the plunge can be viewed as that of a test particle moving in the 
fixed Schwarzshild geometry of mass m. The contribution from the plunge 
phase is estimated using the PN formulae assuming they are valid even 
beyond ISCO. Since the PN representation is usually not reliable inside 
ISCO, this should be a source of error and in general this computation is 
only a crude estimate. Further, the the recoil velocity accumulated during 
the two phases (inspiral and the plunge) can be obtained by taking a vector 
sum of the two estimates. This is achieved by matching the circular orbit 
at the ISCO to a suitable plunge orbit. The final results of the numerical 
estimates of the recoil velocity is presented in Fig [l] We compare the 
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Figure 1. Kick velocity imparted to the remnant of the compact binary 
coalescence as a function of symmetric mass ratio. This figure incorporates 
the contribution from the plunge phase of the binary evolution. The results 
of 2.5PN+plunge is compared against that of Ref. [5] where the estimation 
was done with a 2PN+plunge model. 



recoil velocity based on our 2.5PN inspiral + plunge model with that of 
2PN inspiral + plunge model of Ref. [5j. (Specifications are same as those 
of Fig. 1 of [5].) As is evident, our 2.5PN inspiral + plunge model predicts 
smaller recoil velocity for almost all values of v. The maximum of the curve 
drops from ^ 243 km/s of [Sj to ~ 180 km/s. This may be attributed to 
the oscillatory convergence of the PN series observed in various contexts 
(see e.g. [9]). While our estimates are within the error bars given by [5], 
one should keep in mind that there are also error bars associated with 
our estimates due to systematic effects such as neglect of higher order PN 
contributions. 



5. Conclusion 

Using the MPM formalism, we computed the linear momentum flux due to 
GW emission from inspiralling compact binaries moving in circular orbit 
up to 2.5PN order. This complements the computations of energy flux [2], 
waveform and polarizations [JIHj at 2.5PN order for circular orbits. Using 
the PN linear momentum flux and an analytical model of plunge [5], we 
have estimated the contributions to the recoil from the inspiral and plunge 
phase, as a function of symmetric mass ratio. 
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